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Abstract 

A model of randomly distributed overlapping spheres of different radii is represented to 
describe a heterogeneous porous medium. Two-particle correlation function of the relative 
position of pores of different radii in the medium space was calculated and detailed analysis 
was carried out. The model allows to characterize the disordered porous medium with the 
number of nearest neighbors, the area of all mouths that connect pore with the neighboring 
pores. These paramaters depend on the pore size and porosity and — in addition to the 
specific surface area, porosity and percolation threshold and the distribution function of 
pore size. 
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Introduction 



To describe the dispersion of fluid and interaction of fluid nanoclusters in a disordered 
porous medium is necessary to use a certain model of a porous medium. A detailed description 
of the structural features of the porous body is very difficult and time-consuming. This complete 
information, even if it was obtained is unnecessary. To describe the processes some homogenized 
structure parameters are enough. The fact that linear dimensions of the porous body, with rare 
'"O ! exceptions, greater than the average pore size. In this case, to describe the different properties 
^ \ of the porous body are important different characteristics of the geometry of the pore space. 

In the literature, as a disordered porous medium model consisting of distributed in a certain 
way in the space of inclusions, such as circles, ellipses, cylinders, spheres, ellipsoids [1] is often 
used. These inclusions can either overlap or not. The size distribution can be both monodisperse 
and polydisperse [2]. The most common inclusions are spheres [3J. In the spheres are not allowed 
to overlap, then model is called a fully-impenetrable or hard spheres. In addition to modelling 
of a packing of particle, the model has been used in the study of variety phenomena, such 
as powders, cell membranes jl], thin films [5], particulate composites, colloidal dispersions [6]. 
O The model becomes quite general if spheres are allowed to overlap [7]. The intersection of 
^ \ spheres does not reflect the true physical entity, but can be a way to obtain a disordered 
\ heterogeneous porous media. The distribution of spheres is characterized by a large number 
of different fundamental microstructural properties, including n-point correlation function and 
■ the size distribution function [HIE] - 

In this paper, well-known model of randomly arranged overlapping spheres [10] is general- 
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ized for the case of pores (spheres) of different size, provided a narrow pore size distribution 
(AR/R < 3, where AR — half-width of the pore size distribution, and R — mean radius). We 
calculate the two-particle correlation function mutual arrangement of spheres of different sizes 
in the space environment. In this model, a disordered porous medium can be characterized 
depending on the pore size and porosity of the new options — the number of nearest neighbors, 
the area of all mouths that connect pore with neighboring pores. These options are in addition 
to the specific surface area, porosity, and percolation threshold of the distribution function of 
pore size. With the help of additional parameters is possible to describe the thermal effects [TT] 
and phenomena of nonoutflow |12| in the system consisting of a porous medium immersed in a 
non-wetting liquid [ 13 J . 
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2 A formal approach to the description of two-phase sys- 
tems 



Consider a formal approach to the two-phase system and a porous medium is the two-phase 
system. This approach does not rely on any structural idealization and is applicable to any 
structure of porous bodies. 

Let us designate the space of pores (voids) by m pore and space of solids as m so ^- Then 
total space of the porous material can be written as m = m pore + m so ^- Let us consider the 
function g(x), which fully characterizes the geometry of the two-phase medium [Tj 
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Here, vecx is the radius vector of the point of the space m. Analytical form of function g{x) 
is possible only for materials with a regular structure, such as regular packing of hard spheres. 

Measurable structural characteristics of the medium can be obtained by averaging of various 
quantities containing the characteristic function g(x). Averaging an arbitrary function A(x), 
defined in the space m, means integration over the m space with subsequent divided by its 
volume V. The obtained average value of (A(x))v is equal to 

(A) v = \- f A(x)d 3 x 



V _ 

The average value of the characteristic function g{x) is the ratio of the pore volume to the 
total volume, i.e. porosity <p: 

v = (g(x))v = 7f / g(x)d 3 x (l) 

The characteristic function allows to introduce such characteristics for the porous medium as 
the number of nearest neighbors and the area of all mouths that connect a pore with neighboring 
pores. Such information about the geometry of pores space is contained in the two-point 
correlation function g^O^i)- It can be obtained by the following averaging |15) : 

0i (si) = (g(x)g(x + xx)) v (2) 

The two-point correlation function ^2(^1) is the probability that both ends of the vector x±, 
accidentally thrown into the space m, lie in pores space m pore . By the expression ([T]) it can be 
shown that the function gi(x\) satisfies natural limit relations: 

9 ^ = {^, \x,\ ^00 (3) 

At \xi\ — > 00 the correlation function (j5J) is equal to the product of two independent char- 
acteristic functions with no correlation between them. Therefore the expression (j2J) can be 
replaced by the following expression: 

= (g(x)) v (9(xi)) v , (4) 

Using flTJ one can write the expression (j3J) as the square of the porosity, because averaging 
occurs in the same volume V and each of the factors is equal to the porosity phi. When 
|xi| —> the expression (J2]) can be written as ^ 2 (^?i->o = (fi'(^))y; that, according to ([T]), gives 
porosity <fi. 



2 



3 The correlation function of the relative position of equal 
sizes pores 

Let us consider a medium consisting of randomly overlapping spherical voids of concen- 
tration n. This means that voids — spheres which are randomly placed in space without any 
cross-correlation. The probability of finding the center of a randomly selected empty space at 
a certain point does not depend on location of other centers. The intersection areas are the 
throats where menisci are formed as the liquid fills the pores. This model is called the model of 
randomly placed spheres (RPS). The simplest model operates with spheres of the same radius 
(monodisperse distribution). For the first time such model was suggested by H. Weisberg [7] 
to describe the dispersion of porous media and by W. Haller [16] to describe the structure of 
porous glasses. 

Since a real porous media always has a certain pore size distribution, in the simplest case an 
average radius R and variance DeltaR, we will consider spheres of different radius within this 
RPS model. The introduction of the pore size distribution leads to the necessity to analyze the 
correlations in the mutual arrangement of pores with different radii. In particular, correlations 
are that a large pore has more nearest neighbors than a small pore. As it is shown below, 
in contrast to the situation with the spheres of the same radius, the correlations significantly 
change the connectivity of pores with the neighboring ones and lead to the dependence of the 
nearest neighbors and the surface area of meniscus in pores throats on the pore radius R. 

In the context of the porous medium model under discussion the porosity varphi is the 
ratio of spheres volume to the total volume of the medium. It is obvious that the porosity is 
the probability that an arbitrarily selected point in the space environment is situated in the 
pore medium. In the RPS model this probability equals the probability that there is at least 
one center in a sphere of a certain radius R. 

Let us find the probability P(Vq) that the volume Vo(Vo < V) does not contain any centers. 
This means that all nV centers are located in the volume V — Vq (where n is the concentration 
of spherical particles, V is the total volume of the medium). For an arbitrarily chosen center 
outside of Vq, the probability equals 



Then the probability that the volume Vq contains at least one center can be written as 



As we consider the case when the total number of particles nV is large and the product uVq 
is finite, using the definition of e as the limit of the expression (1 — a) -1 / a for a toO, one can 
write this probability as 



From the equation (j5J) at Vq = 47ri? 3 /3 (the average volume of one sphere) we can get the 
value of the porosity as 



Using the relation (j3J) one can obtain an explicit form of the correlation function, which 
allows to calculate the number of nearest neighbors, the area of meniscus and, accordingly, the 
pore connectivity coefficient with surrounding pores. 
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Let us calculate the pair correlation function g 2 in a system of randomly placed spherical 
pores radius R and R\. This situation is shown in Figured) the considered pore of the radius 
R intersects with pore radius R\. Denote the distance between centers as k and the respective 
heights of the segments formed by the intersection as x and y. Then the correlation function 
g2(R,Ri,K) is equal to the probability that at least one center is located inside the locus of 
points, disposed on a distance less than R (and/or R±) from one end of the segment k in length. 
If k is greater than R\ + R, this area consists of two spheres of radii R and R± and its volume 
is equal to (4/3)-7r(i? 3 + Rf), respectively. If k is less than Ri + R, this area represents two 
intersecting spheres of radiuses R and the distance between the centers of which is equal 
to k. In this case using simple geometric representations one can see that the area volume is 




Figure 1: Schematic illustration of the intersecting pores of different radii R and Ri 

Using these expressions ([2]),©, we can write the correlation function g 2 (R,Ri,n) in the 
generalized RPS model for spheres of different radii as 

( ip 2 k> (R + Ri) 

gi {R,R U K) = < ^ R s +R 3_32 {Ri _ §) _ 32iR _^ (7) 
I if i k < (R + Ri) 

here x = R ~^~ Rl ^ an d y = R + Ri — x — k = Rl ^ — are determined by the geometry 
of the intersection of two spheres (see Fig. 1), k is the distance between centers. 



4 The number of nearest neighbors 

In this case the number of nearest neighbors of pore of radius R in the randomly placed 
spheres model is calculated as the first coordination sphere integral of the correlation function 
([7|). With a glance to normalization we can get 

I r\R+Ri\ 

z{R,R 1 ) = —— d Z Kg 2 (R,R 1} K), (8) 

¥ V por JlR-R^ 

here V por is the volume of one neighboring pore of radius Ri (see Fig. [TJ 
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In the framework of the mean field approximation for the particular porous medium by 
using a function of pore size distribution f(R) we can make the following averaging: 



m = / fm-w- / 

JO VVpor 




Kg2(R, Ri, n)dR 
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Let us consider the number of nearest neighbors for two fundamentally different cases: when 
a small pore is surrounded by other small pores and when a large pore is surrounded by the 
pores of arbitrary size. Moreover, small and large pores will be selected according to the pore 
size distribution. 

Denoting R s and Rb as the radii of the small and the large pores, respectively, from the 
expression (jHJ)we can obtain the number of nearest neighbors for the large pore surrounded by 
pores of arbitrary radius: 



and the number of nearest neighbors for the small pore surrounded by pores of arbitrary 



In figure [2]the dependences of the number of nearest neighbors for the large and small pores, 
surrounded by identical pores of radius R, on this radius R are depicted. From the initial 
assumptions of the model and the calculations it follows that the number of neighbors for the 
small pore surrounded by pores of the same small radius and for the large pore surrounded by 
pores of the same large radius coincides and corresponds with the number of nearest neighbors 
in the identical radius pores approximation. Moreover, as it was expected, one can see that 
for the small pore the number of nearest neighbors falls as the radius of surrounding pores 
increases. For the case when the neighboring pores of considerably outmeasure the considered 
pore, the number of neighbors tends to zero. If the considered pore is large the number of 
nearest neighbors increases as the radius of surrounding pores grows. 

With a glance to (J7J) the expression (jSJ) and therefore (jUj) and ([ID]) depend on porosity. In 
Figure [3] the dependence of the number of the nearest neighbors, calculated according to the 
formulas and ( I10p for different values of porosity is shown. One can see that as the porosity 
increases, the number of nearest neighbors also grows. So in the case of identical pores at 
porosity near the percolation threshold (p ~ 0.2 the number of the nearest neighbors equals 
z ~ 0.2. At cp — 0.4 z = 3.6, at ip = 0.6 z = 5.1 and at tp — 0.8 z = 6.6. When porosity 
tends to one the number of neighbors for the identical pores tends to 8. At first glance it is a 
rather strange result given the fact that in case of the dense packing of the balls, the number 
of neighbors is 12. This difference is due to the condition of an open contact of pores, which is 
ensured by the condition of intersection rather than by the point areas contact. 





radius: 




(10) 



5 




Figure 2: The dependence of the number of the nearest neighbors on the pore radius for 
two pores with different radius: R s (dashed line) and (solid line). Curves are plotted in 
accordance with the formulas d7J), fl5J), (ITUj) . while the porosity is ip = 0.4, R s = 2 nm, and Rf, 
= 4 nm. Straight dotted line indicates the number of the nearest neighbors in the model of 
the spheres of same radius. The intersections with the curves correspond to a situation where 
the pores of R s and Rf, radii are surrounded by the neighbors of the same radii R s and Rb, 
respectively 
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Fisrure 3: The deDendence of the number of the nearest neighbors on the Dore radius for two 



5 Mouth pores area 



In the model of intersection of two spheres, when a pore with radius R is filled, area of the 
meniscus (a pore throat) will be defined as the surface area of the spherical segment, formed 
by the intersection of two pores. From simple geometric considerations area of the segment is 
written as s m i(R, R±, k) = 2nRy (see Fig. [1]). In this case, after averaging over the distance 
between pores k, the meniscus area can be written as 



s m [R, Ri 



1 r R+R ^irR(R{ 



IK 



Rf) 



■ Ann • d,K, 



(11) 



here V = f {(R + R{f — \R — Ri\ 3 ). 



In the framework of the mean field approximation for the particular porous medium by 
using the function of pore size distribution f(R) we can make the following averaging: 



s m (R) 



1 rR+RinRfRl 
f(Ri)dR 1T7 1 



V 



Rf) 



Attk 2 ■ dn, 



H 



(12) 



'0 v J\R-Ri\ 

In Figure H] the qualitative dependence calculated by ()12|) is shown. As the radius of the 
considered pore grows, the value reaches a constant value which is less than a half of the filled 
pore surface area. 
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Figure 4: The dependence of the meniscus surface area s m (R) on the considered pore radius. 
The dependence is calculated by (IT2|) for two overlapping pores in mean-field approximation. 
For the function of the pore size distribution we take Gaussian normal distribution model with 
the average radius of 4 nm and the variance of 0.7. The porosity is tp = 0.4 

Similarly to the dependence of the number the of neighbor pores let us consider the de- 
pendence of the meniscus area of the for a small and a large pores intersecting with a pore of 
arbitrary radius. The respective expressions for the meniscus surface area can be written as 



l Trfl (Rl - (« - Rif) 

Sm\ttb, ril) — ~ 
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(13) 
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Since the meniscus area in the model is defined as the surface area of a spherical segment, if 
the filled pore of radius R intersects the pore of radius R\, the area of the meniscus is not the 
same as if the filled pore of radius R± intersects the pore of radius R. In Figure one can see 
the dependences of the meniscus surface area on the radius of the pores, which they intersect, 
for two different pores. Dependences were calculated by the expressions (|T4|) and ( JT3|) . As the 
radius of the considered pore grows, the value reaches a constant value which corresponds to 
a half of the filled pore surface area. However, such cases are almost impossible, because for 
Gaussian pore distribution they could be observed only when a small filled pore intersects with 
an empty large pore, but the probability of such an event is extremely low and is estimated as 
4 %. 
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Figure 5: The dependence of the meniscus surface area on the radius of a neighbor pore for two 
different pores of radii R s (dashed line) and Rb (solid line). Curves are plotted in accordance 
with the formulas (fT3|) . (lT4"j) . The porosity is (p — 0.4 (g), R s = 2nm Rb = 4nm 

It should be noted that the expression (1TT1) . and therefore (TT2|) . (TT3|) and f[T4"|) do not depend 
on the specific number of spheres n and therefore porosity <p. This, at first glance, strange cir- 
cumstance is connected with the specific character of the model used, namely, with uncorrelated 
random distribution of voids centers in the framework environment. 

Another characteristic of the porous medium is the connectivity coefficient 77. The connec- 
tivity coefficient for a pore is calculated as the ratio of the all menisci surface area to the pore 
surface area. In accordance to its physical meaning it can not be greater than one. With a 
glance to (|S]) and ( II ip . the connectivity coefficient can be represented as: 

1 [ R+R i 7rR(Ri-(K-R) 2 ) , 2 
rj{R, R x ) = —— / LL . g 2 ( R: Ru k )4ttk 2 ■ dn (15) 

bVVporP J\R- Rl \ K 

here S = 4irR 2 is the total surface area of the filled pore. The expression (jl5[) depends on 
the radius of the considered filled pore and the radius of a neighboring pore. In the framework 
of the mean field approximation we make an averaging over the neighboring pores using the 
function of the pore size distribution f(R): 

V(R) = —— / / LL . f{R l )g 2 {R^ R u k)Attk 2 ■ dndR 1 (16) 

VVp 0r (p Jo J\R- Rl \ K 
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Numerical analysis of the connectivity coefficient expression ( 116]) by the correlation function 
(J7J has shown that the dependence of averaged connectivity coefficient over neighbors on the 
radius can be written as: 



V (R) = Vo /R a , 



(17) 



here 7]q is some constant and a is an exponent lying in the range from zero to two. For 
a porous medium Libersorb 23 with porosity of 0.56 and Gauss pore size distribution: mean 
radius R = 4.2 nm and dispersion 0.6 nm in Fig. [6] is shown the dependence of connectivity 
coefficient rj on the radius of the pore R. The exponent a is 0.4. 
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Figure 6: The dependence of the dependence of connectivity coefficient r) on the pore radius in 
the mean field approximation for the porous medium medium Libersorb) 



6 Approximation of the monodisperse distribution model 

Let us show that when the pore size distribution is a delta function, i.e., actually all randomly 
placed spheres have the same radius, the obtained expressions are the ones known from jTUJUT]. 
Indeed, if we replace Ri by R, the two-point correlation function g\ is in fact equal to the 
probability that at least one center is located inside the locus of points disposed on a distance 
less than R from one end of the segment k in length. If k is greater than 2R, this area consists 
of two spheres of radius R and its volume is equal to (8/3)7ri? 3 . If k is less than 2R, this area 
consists of two intersecting spheres of radius R, the distance between the centers of which is 
equal to k. When R = Ri, from the equations (jTJ) we can get 



x 



R 2 -(k- Rf -k 2 + 2kR 
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2k 



2k 



+ R 



y = 2R-x-K = 2R+--R-K = R-- 
y 2 2 

Therefore, in ([7]) x — y and the volume of the intersection area of identical spheres is 

3k k 3 \ 4 , . . 

1 + —ttR 18 

AR 16R 3 J 3 K ' 

Substituting this expression (I18p in the correlation function ([7]) we get: 
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g 2 {R, k) 
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and finally: 



k<2R 



1 J_ 3k K 



k < 2R 



(19) 



(20) 



The resulting expression (1201) satisfies the limit relation ([3]) and coincides with the expression 
for the correlation function from jlOj . 

The number of nearest neighbors in the framework of the equal radius spheres model is 
calculated as the integral of the correlation function in the first coordination sphere. With a 
glance to the normalization of the correlation function we can obtain: 



z(4>) 



2R 



(21) 
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The resulting expression of the number of nearest neighbors only depends on the medium 
porosity and is independent from the pore radius, which is clear from the fact that in this model 
the radius variation actually leads only to the variation of the system scale. The number of the 
nearest neighbors dependence on the porosity (j21j) is shown in Figure [71 As it was expected, as 
the porosity increases, the number of neighbors grows. Note that the RPS model predicts the 
percolation threshold, which corresponds to the formation of a infinite connected cluster. The 
formation of an infinite cluster at an average number of throats per one pore equal to z ~ 2 is 
possible if the porosity is close to the percolation threshold (exceeds it) tp c = 0.29. And it is o 
shown in [18] that at such z ~ 2 the part of pores available for filling in a disordered porous 
medium tends to unity. 




Figure 7: The dependence of number of the nearest neighbors on the porosity tp for the RPS 
model with spheres of the same radius 

Considering the meniscus as a spherical sector, the meniscus area can be written in the 
following form 



s m (R, k) = 2kR(R-k/2) 



(22) 



Consideration of the specific meniscus form leads to the appearance of a numerical factor 
in this expression and does not affect the character of the obtained dependencies. 
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Then the connection coefficient, which is the ratio of the area of all menisci to the surface 
area of the pore: 



1 
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4nR 2 V (p 



s m (R, K)g 2 {K) 



(23) 
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Figure 8: The dependence of the connectivity coefficient r\ on the porosity (p for the RPS model 
with spheres of the same radius 

One can easily see that the expression ( !2"3"j) is independent from the pore radius, and depends 
only on the porosity of the medium. The dependence of rj on the porosity is shown in Fig. |HJ 
Note that when the porosity is </? — > 1, the ratio of menisci surface area to the pore surface area 
must tend to unity 77—7-1, which is observed in Fig. [HJ 

In the RPS model in case of identical radius pores the number of nearest neighbors can also 
be calculated in the next approximation. The structure of the pores in the porous medium can 
be represented as a random lattice whose nodes are voids and links are throat. By the number 
of nearest neighbors we mean the number of links radiating from a node (the number of void's 
throats). Let us find the void distribution function over the number of the nearest neighbors. 
Formally, the model assumes the existence of a void with any coordination number, but too 
small and too large values of z are improbable. The probability p z that a randomly selected 
void has strongly z throats is equal to the product of the probability that out of (n — 1) voids 
centers, located in the unit volume containing the selected void, exactly z are situated in a 
sphere of radius 2R with the center in the selected void and the probability that the remaining 
(n — 1 — z) voids are situated outside the sphere of radius 2R with the center in the selected 
void. Combinatorial factor allows taking into account the variants of the arrangement of n — 1 
throats by the number of the nearest pore neighbors: 




(Sa) z exp-Sa/z\ 



(24) 



The average number of the nearest neighbors z is equal to 




(25) 



The calculation of the amount with a glance to fT24|) gives 
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The comparison of the number of the nearest neighbors computed by the coordination 
function and by the probabilistic approach in the random lattice model is shown in Fig. |9j 
For the porosity up to <£>~0.5 two dependencies are close, but for the high-porous medium an 
average number of the nearest neighbors computed by (|25j) significantly increases resulting in 
nonphysical behavior of the connectivity coefficients rj, which becomes greater than one for 
the porous media with porosity <fi > 0.8. Use of this approach in case of low porosity can 
be justified by analytical expressions for the connectivity coefficient rj and the number of the 
nearest neighbors. 
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Figure 9: The dependence of the number of the nearest neighbors on the porosity <p for the 
RPC model with two spheres of the same radius calculated by the formula f[2"Tj) (solid line) and 
by the formula ()25p (dashed line) for the the random lattice model 

However, the model consisting of the same radius spheres as well as the random lattice 
model lead to such a behavior the " nanoporous environment - non-wetting liquid" system 
when the fluid either fully outflow or does not outflow at all depending on the filling degree 
and the surface energy. It contradicts the observed experimental data [12J, in which one can 
observe a partial leakage (not all the liquid flows out of the porous medium). Therefore, for 
description of the inflow and outflow of the liquid from the porous medium the generalized RPS 
model will be used. To describe mechanical energy dissipation and thermal effects we will use 
the same radius RPS model. 
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